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We consider an N-server queuing system with Poisson arrivals and 
exponential service, in which arriving customers must pass through a 
gate into a waiting room before becoming eligible for service. Customers 
who find the gate closed wait outside until the gate opens; customers 
inside the waiting room are served at random. When the last customer 
inside acquires a server, the gate admits all those outside and then closes 
again. If no customer is waiting outside when the gate opens, the gate 
remains open until there is a queue of k waiting customers. 

Service offered by this system is intermediary between random service 
and order-of -arrival service. As long as the gate is open and fewer than 
N + k customers are in the system, service is purely random. The param- 
eter k can be regarded as a threshold at which the queue is judged too long 
to permit random service to continue. 

Our main results are (i) the Laplace-Stieltjes transform of the equilib- 
rium distribution of the waiting time of an arbitrary customer and (ii) 
a comparison of the second moments of the waiting time for different 
values of k with those of the waiting time under random service and order- 
of -arrival service. The service is shown to be "nearly random" at low 
loads and "not quite order-of-arrival" at high loads; for higher values of 
k this transition occurs at higher traffic intensities. 

I. INTRODUCTION 

Switching systems, particularly electromechanical switching systems, 
are often constructed so that if several customers are awaiting service 
simultaneously, they will receive service in what is essentially random 
order, i.e., a server which becomes idle will choose its next customer 
at random from the queue of customers awaiting service. Such an 
arrangement may be satisfactory when the traffic intensity is low. 
but as the intensity increases, a progressively greater number of 
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customers will have to wait an undesirably long time, and the quality 
of service may become unacceptable. 

There are available, however, methods of providing service other 
than "random service" ; the most obvious one is service in order of 
arrival. The quality of service, which depends in part on the variance 
of the waiting time, will still diminish as the traffic intensity increases, 
but not as quickly as when random service is used. In fact, order-of- 
arrival service is the "best" discipline (at least when the order of 
service does not depend on individual service times) in the sense that, 
for a given traffic intensity, and hence mean waiting time, the variance 
of the waiting time is smallest. 1 Unfortunately, it may not be worth- 
while (or even possible) to build a system which offers service in order 
of arrival. One is therefore led to consider an intermediary queue 
discipline, one for which the variance of the equilibrium waiting time 
lies between that of random service and that of order-of-arrival 
service. 

Suppose we have an M/M/JV queuing system with customers arriv- 
ing at rate X and requiring a mean service time n~K Suppose further 
that customers must first pass through a gate into a "corral," or 
"waiting room," before becoming eligible for service. So long as there 
are not more than N customers in the system, the gate remains open ; 
an arriving customer enters the corral immediately and, if some server 
is idle, begins service. As soon as a customer has to wait (having 
found N customers in the system), the gate closes until that customer 
enters service. The gate then opens to admit all those who have 
accumulated outside the corral, and immediately closes again, re- 
maining closed until all those inside have acquired a server, and so on. 
Thus the customers are admitted in bunches to the corral, and once 
inside, they are served at random. It is assumed that the corral has 
an unlimited capacity. If there is no customer waiting when the gate 
opens, the gate remains open until there is again a customer who has 
to wait. 

A gated queuing system merits consideration, not only because of 
the anticipated improvement over random service, but also because 
it can overcome design problems in certain telephone equipment which 
might otherwise lead to poor service for some customers. For example, 
in some switching equipment, each server hunts in a fixed sequence 
over a group of customer-sources; when a customer is found, the 
server stops to provide service and then resumes hunting from that 
point. Such a hunting procedure may be desirable from a hardware 
viewpoint, but it can result in uneaual service among customers. 
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Gates have been successfully used to improve service in the manner 
described above, that is, by temporarily blocking subsequent bids 
for service until all those customers present have been served. 

Some discussion of the gated queuing sj^stem is given in a 1953 
paper by Wilkinson. 2 Theoretical results were summarized in 1958 
by Beckwith, 3 although he gives few details as to how the results are 
derived. The model we consider in this paper is more general than the 
one described above. We shall assume that as soon as there are k 
customers (rather than one) waiting, the gate closes until all k cus- 
tomers have entered service. Thus the parameter k can be thought of 
as a threshold at which the queue is judged to be too long to permit 
purely random service to continue; the system enters the "gating 
mode," and the gate admits customers in bunches as described pre- 
viously. If the gate opens and there is no customer waiting outside, 
then the system leaves the gating mode and the gate remains open 
until there is again a queue of k waiting customers. 

In Section II we obtain a recurrence relation for the probability- 
generating function of the nth bunch-size and, after that, the generat- 
ing function and moments of the equilibrium bunch-size distribution. 
Using these results we determine in Section III the Laplace-Stieltjes 
transform of the distribution of the equilibrium waiting time of an 
arbitrary customer. In Section IV we obtain the first two moments of 
the equilibrium waiting time, and we make some comparisons of the 
second moments of gated service for different values of the parameter 
k with those of random service and order-of-arrival service. 

In several places in the text, we specialize a general result to the 
case k = 1, since that is the simplest of our gated systems. This 
allows us to verify that our results then agree with Beckwith's, and it 
often reduces a complicated expression to one that is more easily 
comprehended. 

We can assume, without loss of generality, that the system is 
initially empty. Since we assume that the traffic intensity, p = X/iV/Zj 
is less than unity, the number of customers will always return to zero 
in a finite length of time, regardless of how many customers may be 
present at the beginning. Consequently, our equilibrium results do 
not depend on the initial conditions. 

Several other authors have considered systems which operate in a 
manner similar to ours, but they all take k = 1 and assume that the 
underlying system is an M/G/l queue. The "generations" in the 
M/G/l queue, as described by Kendall 4 and later studied by Neuts, 5 
correspond to the bunches in our model. Nair and Neuts 6 ' 7 subse- 
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quently consider the waiting time distribution for the M/G/l queue 
under the assumption that the queue discipline was either longest- 
processing-time-first or shortest-processing-time-first. 

II. DISTRIBUTION OF THE EQUILIBRIUM BUNCH-SIZE 

Let X„ be the number of customers in the nth bunch; that is, there 
are N + X n customers in the system the instant after the gate closes 
for the nth time. Thus X n is the number of customers waiting outside 
if that number is positive, and is k if there was no one waiting, when 
the gate was opened for the (n — l)st time. Let T n be the time it 
takes to serve X n customers. Then, denoting by p d (T n = t) the 
density of T n at t, we can express the distribution of X„ in terms of 
that of X„_i by 

P(X n = j) = £ P(X n -i = i) / " p (*« = il X »-i - *'» T - 1 = 

■p*(ZV, = *|X„_i = i)dt. (1) 

But, given that X„_i = i, Tn-i has a gamma distribution with param- 
eters i and Nn ; and 

P(X n = j|X„_i = i, T n _! = t) = P(X n = j\T n -l = t) 

M-%-^ + 6-^, j = k. 

The extra term e~ xt in this expression is the probability that no one is 
waiting when the gate opens ; when this happens, the next bunch-size 
is necessarily k. Substituting these expressions into eq. (1) gives 

W, = i) = |p(x.-, = i)( r ^)' 

x[(r^)'( i+ r x ) + 4 

where p = \/Nfx and 8 jk is the Kronecker delta. Now, letting 
Pn(u) = ±P(X n = j)u' 

7-1 

be the probability-generating function of the distribution of X n , we 
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have 

'.(•) - r~* ( i + pa-,) ) + <«* - W~ (ih} (2) 

Starting with P\(u) = u k , we can determine successively the P n (u). 
When k = 1, eq. (2) agrees with Beckwith's eq. (1). 

We wish to obtain the distribution of the equilibrium bunch-size 
X = limn-oo X n . To see that an equilibrium distribution exists, we 
observe that the bunch-sizes X n form a Markov chain which is ir- 
reducible and aperiodic. Since we are assuming p < 1, the number of 
customers in the system cannot grow without bound, so the states are 
not all transient or null. Therefore all states are ergodic, and there is 
a unique equilibrium distribution. 8 The distribution of X = lim„-.oo X n 
has a probability-generating function P(u) = lim„^. x P n (u), which 
satisfies 

f(»)-f( 1 + < , ( 1 1 - u) ) + («'-i)f(rT7> 

This can be written in the form 

P[r(w)] - P(u) = F(u) 
by setting 

r(u) = 1/[1 + P (l - w)] 
and 

F(u) = (1 - u*)P[l/(l + P )]. 

The solution to this functional equation is 9 

P(u) =v~ t Flr n (u)l 

n=0 

where r„ is the nth iterate of r and -q is a constant. To evaluate this 
expression, we first need to find r n (u). We have r (u) = u and 

- = J + /_ ^ (3) 

Letting y n be such that 

rn = Kst! + 1 + l 

Vn P 

transforms eq. (3) into a linear homogeneous difference equation, 
which is easily solved. Thus we determine that 

r (u) = 1 ~ PU + (M ~ 1)p " 
nW 1 - P u + (u - l)p»+ l> 
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The solution to our functional equation is therefore 

Since P(l) = 1, we must have y = 1 ; since no bunch-size can be zero, 
P(0) = 0. This determines P[l/(1 + p)] f so that finally 

v ji _ r i - /»« + (« - i)p n t 

FW — ih^ti 

~ x MO)' ( 4 ) 

where /i(w) is the numerator in the second term of P(w) above. By 
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Fig. 1 — Possible states of the system and the transition rates between states. 
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setting k = 1, we can reduce eq. (4) to Beckwith's result, his eq. (2) : 



P(u) = 



1 

0(1) 



where 



u — 

1 - 


-;♦•<»-(.::: 


\ —. •> 


. £ p b 



n f 1 ~ P^hl 



The moments of X are obtained by differentiating P(u). We find 
that 



and 



Var (X) = 



ft 



fc p 



ft 



M0)(1 -p) .1 +p 1 -P A(0)(1 



=s] 



III. DISTRIBUTION OF THE EQUILIBRIUM WAITING TIME 

Let W be the waiting time to the point of entering service of an 
arbitrary customer when the system is in equilibrium. The distribution 
of W depends on which state the system is in when the customer 
arrives; i.e., it depends on the number of customers already in the 
system and on whether the gate is open or closed. These states, to- 
gether with the transition rates from one state to another, have been 
enumerated in Fig. 1. Let 

Pj = P[_j customers in the system], j ^ 0, 

p c j = P[j customers in the system, gate closed], 

j = N 4- 1, N + 2, • • • , N + k - 1, 

P°j = PU customers in the system, gate open], 

j - N + 1, N + 2, • • • , N + k - 1. 

Obviously, p) + p° = Pj- It is also clear that when j ^ N, the gate is 
open, and that when j 1= N -\- k, the gate is closed. The values of the 
Pi are just the equilibrium state probabilities of an M/M/N queue, 
which are 



Pi = ^-jr- Po = -jr po, 

A P N) 

PN+j = PPN+j-l = P J PN - P J — jyj 



j - 0, 1, • • • , N 
3>0 



Po, 
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To find p% we equate the rates at which the system leaves and enters 
the state [j customers in the system, gate open}. From Fig. 1, we see 
that 

(1 + p)P°n+j = PPn+j-i + pSr+i+i, i = 1, 2, • • • , fc - 1, 

with p° N+k = and p° N = p N) where p N is known from the above. 
The solution to this equation is 

Px+i m Y^JkVN, j = 0, 1, • • • , k - 1. 

It then follows that 

Pn+j = PN+i - Pn+j 

- rrr^VpAT, j- l, 2, ..-,/c-i. 

To find the distribution of the waiting time W, we shall consider 
what happens when an arriving customer encounters one of the 
following conditions : 

Hi : < N customers in the system, 
Hi : the gate is closed, 
H%.j : N + j customers in the system, and the gate is open, 

j = 0, 1, • • • , k - 1. 

From the above computations, 

P(Hi) = 1 - £ p j Pn = 1 - t-=-Pn> (5) 



y=o 



1 - P 



P(Ht) = *£; p% +} + £ p™ = t-^tt* p* ( 6 ) 

y=i ;=* i p 

and 

*>(#•.*) = YlA V* 3 = 0, 1, • • • , k - 1. (7) 

An arriving customer who encounters Hi immediately gains access 
to a server, so 
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Fig. 2 — Relations among the variables used. 



An arriving customer who encounters H 2 will have to spend some 
time, Z, outside the gate waiting for the gate to open (see Fig. 2) ; 
once he gets inside, he will then have to spend some additional time, 
Y, waiting to be chosen for service from the bunch he entered with. 
The total amount of time the customer spends waiting to enter service 
is therefore W = Z + Y, where Z can be regarded as the residual 
lifetime of an interval Go during which the gate is closed. Thus we 
can write 

V d iW = t\H 2 ) = f f X p"(Z +Y = t,G = y,Z = x)dydx 
= r f P d (Y = t-x\G = y,Z = x) 

Jx=0 Jy=0 

■p d (Z = x\G = y) V d (G Q = y)dydx. (9) 

We first need the distribution of Go. Disregarding our arbitrary cus- 
tomer temporarily, let G be the equilibrium length of the time in- 
terval during which the gate is closed. If we are given that j customers 
entered the waiting room when the gate last opened, then the gate 
will remain closed for j service times. Since the equilibrium probability 
that j customers entered when the gate last opened is P(X = j), 
we have 

p*(G = y) = ±P(X = j)b*J(y), 

3-1 

where b (y) = Nn exp ( — N/iy) and the asterisk denotes convolution. 
The mean of this distribution is easily found to be 



E{G) = 



E(X) 



k 



AT M Nn(l - p)h(0) 
Let to be the epoch at which the arbitrary customer arrives, and let 
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Go be the length of the G-interval containing t . Then (see Ref. 10) 

= Inm(i - P )h(0) t P(X = j)h*'{y). (10) 
/o j =1 

Next, the distribution of Z, given G = y, is uniform on (0, y), so 

fl 

p d (Z = x\G = y) = - 



^ x £ « 
»' (11) 



0, .x>y. 
We also have 

p d (7 = « - x|G = y, Z = x) = p d (Y = t - x\G = y) 

CO 

= X) p rf (7 — £ — x\n arrivals in (0, y), G = y) 

■p d (n — 1 other arrivals in (0, y)\G = y) 

where the first factor in each summand reflects the fact that the ar- 
bitrarily chosen customer has probability \/n of waiting one service 
time, 1/n of waiting two service times, • • ■ , and l/n of waiting n 
service times. Substituting eqs. (10), (11), and (12) into eq. (9), we 
obtain 

p d (W = t\H t ) 

= iiMl - p)H0) t £P(X = j) t f b*'{t - x) 

K j=l n = l 1=1 Jx=0 

■ T b*J(y)e-^^p-dydx. (13) 

Notice that the range of integration of the inner integral was reduced 
from (0, oo) to (x, «) because of eq. (11), and that of the outer in- 
tegral was reduced from (0, ■») to (0, t) because b* l (t - x) = for 
x > t. By setting fc - 1, Nn = 1, and X =» p in eq. (13), we can ob- 
tain Beckwith's expression for the corresponding density in his model, 
the density of W given that the arbitrary customer finds more than 
N customers in the system. 

We can calculate the Laplace-Stieltjes transform, \{/(s), of the dis- 
tribution (13), but the algebra is long and tedious. The results are 
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given below, one in terms of P(X = j) and the other explicitly: 
f(Nns) 



(1 - p)M0) 
pks 2 



P(X - n lnor (1 +«)(! +S + Sp) 

P(X - 1} log <l + , + p)(l + ,)-p 



+£jV*-4-(mr-(r#kr 



+ 



( 



1 +s 



(1 +s)(l 4-p + s) 

- r fc P »(i - p) 2 ) / i -p- y- 1 
ALd -p" +1 ) 2 lu -p n+1 / 



^n 



(14) 



1 -P 

/bps 2 

, / (1 + 8 - p - Sp"+')(1 + g - p - Sp"+2) \ 

8 \ (1 - P)[d + S) 2 - p - 5(1+ 8 4- p)p n+1 ] / 

* /fc- i\/ i-p-i \w/ p-i(i- p )« y-i 

Aw-i/\ i-P" / \d -p")(i -P n+1 )/ 

te' 1 f/ 1 -p"+i \» / (l - p«+»)(i + a) \« 

m=l>4\ 1-p / U+S-P-Sp n+1 / 

_ / (l -p"+')(l +S + S p) \"' 

v i + s - P - s P »+ 2 y 

, / (1 - p"+Q[(l + s) 2 4- ap ] \-l 

r V(i4- S )2-p- S (i + s + p) P «+>; j 



(i - P ) - 



s(l + s 4- p)p' 

(1 + s - p - sp n ) k 



I _ p n+l p " (! + s _ p _ sp n +1 )*-l 

1 (1 + g - p - sp" +1 )* 
14-s'(14-s-p- sp"* 2 ) 11 - 1 



+ 



1 [(l 4- s) 2 - P - «(i + s 4- P )p«] fc 



1 4- s [(1 4- s) 2 - p - *(l + s + p)p" +1 ] i - 1 



(15) 



When k is set equal to 1 in eq. (15), the sum from 2 to k is trivially 
zero, and the expression in the last pair of braces collapses to zero; 
the transform then reduces to 

ps 2 n=0 (1 - p n+1 ) 2 



Xlog 



Ld- 



+ 8 — p — sp n+1 )(l 4- s - p - Sp n+2 ) 



p)[(l + s) 2 -p-s(l +S + p )p"+l]J' 

(k = 1). (16) 
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The final portion of the waiting time distribution needed is 
fj(t) = P d (W = t\H 3 .j), j - 0, 1, • • • , fc - 1- Let Wj be the waiting 
time of a customer who arrives to find the gate open and N + j in the 
system. Then Wj has the density //(*), and, letting E\(t) denote the 
exponential density Xe _x ', 

j = 0, 1, ■ • • , k - 2 
/k-i = r -Eatm + r #*£ + " ' " + £ -^m" 

The reasoning behind these equations is that, if an arriving customer 
finds the gate open and j(< k ~ J ) customers waiting, then either 
the next event is an arrival, in which case he waits until that event 
plus an additional time distributed as Wj+i, or else the next event is a 
departure, in which case he waits until that event, after which either 
he is served immediately [with probability 1/0' + 1)] or someone 
else is chosen and our customer must wait an additional time dis- 
tributed as Wj-i [with probability j/(j + 1)]. If, on the other hand, 
the arriving customer finds the gate open and fc — 1 customers waiting, 
then the gate shuts behind him and he waits either 1, 2, • • • , or fc 
service times, each having probability 1/fc. Taking Laplace-Stieltjes 
transforms of the equations, and denoting the transform of /,■(*) by 
<t>i(s), we obtain 

UMNns) = (1 + p + 8)j<f>j-i(Nit8) - (j - l)<f>j-2(NuLs) - 1, 

k j = 1, • • • , fc - 1, (17 ) 

| 0t _ 1 (Ar M8 ) = fc I[i-( r L-)]. 

For any particular fc, this set of equations can be solved explicitly by 
successive substitution. Finally, using eqs. (5), (6), (7), and (8), we 
can represent the Laplace-Stieltjes transform of the distribution of 
the waiting time W as 

4>(Np,s) m 1 - —L_ VN + PN "ff ^^ +/C1M 

L — p j'=0 1 ~ P 

+ irrpp«*Wi»)> (18) 

where <t>j(N^s) is given by eq. (17) and $(Nns) is given by eq. (14) 
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or eq. (15). When k = 1, eq. (18), with the help of eq. (16), can be 
written explicitly as 



<j>(Nns) = 1 
•log 



1 



Vn 



* p , + s^#, r ± 



1 - p r " • 1+ 8 "" ^ S 2 «" „f o (1 - p»+i)» 

;i + s - p - sp n+i )(i + s - P - sp) n ^ 



. (1 - p)[(l + «)« - p - 8(1 + S + p)p 



n+2 T 



(k = 1). 



IV. MOMENTS OF THE EQUILIBRIUM WAITING TIME 

Since the mean waiting time does not depend on the queue dis- 
cipline (see Ref. 11), the mean is the same as for a simple queue with 
service in order of arrival, i.e., 

The second-moment computations are fairly lengthy, finally yielding 
2p N 



E(W 2 ) = 
+ 



JV M (1 - P k ) 
k VN p k ~i{l - p) \(k- \){k - 2) 



!/ !/(p) - p^M'(l) -p*-W(l)] 

2 + 2p + 3p 2 



1 -p 3 



] 



(W(l - P *) I 6 

r2_+3p + V-P 4 l 2-1 4 p + 5p 2 + 2p 3 + p* | . 

^ { } L (1 " P 2 )(l - P«) . (1 " P)d - P 2 )(l " P 3 ) r J 

where Af is a function denned by 

M(x) = - ifzty; (0) = £ x'E(W\H s .i). 

The variance of W is obtained bv subtracting the square of the mean 
of W: 

Var (W) = SOP-) - ( ^f_ p)4 - 

For comparison purposes, we also need the second moments for 
order-of-arrival service and for random service. When service is in 
the order of arrival, we obtain 11 

2p N 



E(W 2 ) = E(W*\W > 0)P(W >0) = 



(20) 



(tfp)»(l - p) 3 
When service is at random, the second moment can be written as 11 



E(W 2 ) = 



VPN 



(AWl - p) 3 2 - p 



(21) 



(22) 
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Observe that the second moments depend on N only through the 
factor p N /(Nfx) 2 , assuming the value of p is fixed. This is true also 
for the second moment in eq. (19), since each of the M-terms con- 
tains a factor (Nn)" 1 . It is therefore convenient to consider the ratio 
of E(W 2 ) for the gated system [eq. (19)] to the second moment for 
the order-of-arrival system [eq. (20)], i.e., 

R{k> p) = MLZ-g)l l p M'( P ) + M( P ) - P *-W(l) - p*-W (1)] 

, fc p *-i(l - P y f (fc- l)(fc-2) [ 2 + 2 P + 3 P 2 1 
+ 2(1 - p k ) J 6 L 1 - P 3 J 

,,. t , f 2 + 3p + 3 P » - p* l , 2 + 4p + 5p 2 + 2p 3 + p* 

+ (A ' " 1} L (i-p 2 )(i-p 3 ) J (i - p)d - P 2 )d " P 3 ) 

Because this ratio is independent of N, it provides a useful tool for 
examining the effect of the value of k on the second moment of the 
waiting time W. Thus we shall be interested in determining its 
properties. 

The function R(k, p) has been plotted in Fig. 3 on the interval 
< p < 1 for a variety of values of k. We observe that R(k, p) is 
bounded from below by unity, increases as k increases, and is bounded 
from above by 2/(2 - p), the ratio of eq. (21) and eq. (20), which 
corresponds to k = ». This general behavior is, of course, just what 
we expected a priori. In order to demonstrate analytically that 

l^R(k )P )^^-, (23) 

I — P 

we first introduce the inequality 
1 j + *^vnvin n< l 1+1 



^E(W\H 3 .j)^ 



^-p-^rr,^^ 2 2 _ p , 

j = 0, 1, ■ ■ • , k - 1. (24) 

The left half of this inequality is demonstrated by considering what 
happens when an arriving customer finds j others waiting, but no 
more arrivals are permitted to enter the sj'stem: the customer's 
expected waiting time will decrease to (j + 2)/2iV M , since the customer 
will have to wait either 1,2, • ■ • , or j + 1 service times, each with 
probability (j + I)" 1 . The right half of the inequality is demonstrated 
by considering what happens when there is no gate at all to block 
future arrivals when a threshold k is reached : the mean waiting time 
E(W | H 3 .j) would then increase to the corresponding mean in a system 
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Fig. 3 — Ratio of second moments. 

employing purely random service. In a random service system, a 
customer who arrives to find j others waiting has an expected delay of 



Nil' 2 



2-p' 



y = o, i, ■■•,*- 1, 



a fact which is derived in the appendix. 

Using the left half of eq. (24), we can substitute (,/ + 2)/2Np for 
E{W\H 3 .j) in R(k, p) and combine terms, obtaining 



R(k, P ) ^ i 



+ v^d - p) 1 2 + 3fc(i - P )(i -p 2 )^ ni - p) 2 i 



12(1 - p k ) 



(1+p)(1+p + p») 1 +P + P 2 



from which it is obvious that R(k, p) =^ 1. Similarly, using the right 
half of eq. (24), we can substitute (j + 2)/iVp(2 — p) in R(k, p). 
Combining terms, we obtain 



R(k,p) ^ 



kp*(l - p) U»(l -p) 2 (2 + 3 P 2 ) 



2 -p (2-p)(l -p*) 12(1 + p + p 2 ) 



, fr(l -p)(2 4-2p + 5p 2 + p 4 ) 2 4- 7p + 10p* 4- 8p 3 4- 3p 4 

"T j/ll W11 I t\ I 



4(1 +p)(l + P 4-p 2 ) 



6(1 + p)(l 4- p + p 2 ) 
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from which it is clear that R(k, p) ^ 2/(2 - p). Thus, eq. (23) is 
established. 

Perhaps the most striking feature of Fig. 3 is that all the curves 
approach the same value, 7/6, as p — > 1. It is easy to demonstrate, 
by using eq. (22), that this must occur. Since the conditional means 
E(W\H 3 .j) remain bounded as p-*l, M(p) -> M{\) and both are 
finite. Thus the first term of R{k, p) approaches zero as p — *• 1. In the 
second term, the factor (1 — p) 4 in front causes all but the last term 
in braces to approach zero. Thus, 

„ n ,x v k P k ' 1 2 + 4p + 5p 2 + 2p 3 + p 4 _ 7 

R{k > 1} - ™1 2(1 +p+ ••• +p^)" (1 + P)d +P + P 2 ) "6 

In order to gain some insight as to why the curves meet at a common 
point at p = 1, we will find it helpful to consider a supplementary 
variable, the fraction of time the system spends in the gating mode. 
When the system is in equilibrium, this is simply the probability that 
an arriving customer finds the gate closed, and is given by eq. (6) : 

kp k 
P (system is in gating mode) = ?-— — - k Vn. 

This quantity has been plotted in Fig. 4 as a function of p, for the 
arbitrarily chosen value N = 7. It can easily be seen (and is intuitively 
obvious) that when p is very close to 1, the system spends almost all 
its time in the gating mode. But when the system is in the gating mode, 
the system's operation is independent of the value of k; it is only 
when the gate is open that the threshold value k can have any effect. 
Thus, as p approaches 1 , the system becomes independent of k ; so we 
can expect the curves in Fig. 3 to be independent of k at p = 1. 

Another feature of the curves in Fig. 3 is that the slope at zero for 
k ^ 2 is the same as the slope of 2/ (2 - p) ; but the slope for k = 1 
is zero. The reason for this becomes clear when we realize that when 
k = 1, order-of-arrival service is guaranteed until there are N + 3 
customers in the system, while fcg2 makes it possible for passing to 
occur as soon as there are JV + 2 customers in the system. For small 
p, N + 3 in the system is much less likely than N + 2. 

The ratio in cq. (22) is convenient because it is independent of N. 
The variance of a distribution, however, is also frequently of interest. 
It is clear from Fig. 3 that the second moments, and therefore the 
variances, increase with increasing k. There is, therefore, a similar 
family of curves, one family for each value of N, obtained by computing 
the ratio of the variance with threshold k to the variance with order-of- 
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Fig. 4 — Equilibrium probability that gate is closed (N = 7). 



arrival service. In Fig. 5 we have plotted some of these curves, again 
for N = 7. It is obvious that these curves have the same general 
shape as those in Fig. 3. The main differences are (i) that the ratio of 
the variances when k = °° (random service) goes to 3 as p — * 1, while 
the ratio of the second moments when A- = » goes to 2, and (it) that 
the ratios of the variances for k < °o go to 8/0 as p — > 1, while the 
ratios of the second moments go to 7/6. These facts are easily verified 
analytically by letting p — * 1 in the actual expressions for the ratios. 
Since the second moments (and variances) increase as a function of 
k, there arises the question of why one might consider a threshold 
value k greater than 1. Clearly, if the queuing systems were otherwise 
equivalent, one would prefer A; = 1 to any larger value of k. But it is 
equally possible that a queuing system would be more costly to operate 
when it is in the gating mode, since more bookkeeping is necessary : 
each waiting customer must be classified as "inside" or "outside" 
the waiting room, and these labels must be changed when the gate 
operates. One would then prefer a system in which the gate were used 
as little as possible (large k), consistent with an acceptable quality 
of service. The resulting tradeoff between cost and quality of service 
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can be resolved only by examination of the particular application at 
hand. 



V. CONCLUDING REMARKS 

In summary, our main result is the specification of the distribution 
of the equilibrium waiting time of an arbitrary customer in a queuing 
system whose service discipline is a compromise between service in 
order of arrival and random service. Our model contains a parameter, 
k, which determines how "close" the discipline is to order-of-arrival 
service or to random service. We have seen that the variance of the 
waiting time is bounded from below by the variance for order-of- 
arrival service, and that as k increases, the variance increases, ap- 
proaching the variance of the waiting time when random service is 
employed. We also found a convenient quantity, R(k, p), which is 
independent of the number of servers, and which, together with Fig. 3, 
allowed us to examine the effect of the threshold k on the waiting time. 
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Figure 3 shows how, for fixed fc > 1, the service changes from "nearly 
random" to "not quite order-of-arrival" with increasing load, and how 
this transition occurs at higher loads as k increases. 

There are, of course, other variables which can be derived from the 
system we have described; for example, in studies of equipment life, 
it might be useful to know the distribution of the number of gate- 
closings that occur in an interval of length t. It did not seem worth- 
while to pursue such questions in the present study, which deals with 
gating from the viewpoint of traffic performance. 

APPENDIX 

Suppose we have an M/M/N queuing system employing random 
service, with arrival rate X and mean holding time p~ l . The mean wait- 
ing time to the point of entering service of an arbitrary customer in 
such a system is 

m = M 

WmU-p) 2 ' 

where p = \/Nn < 1 and pw is the known equilibrium probability 
that there are exactly N customers in the system. We wish to deter- 
mine the mean waiting time, my, of a customer who arrives to find 
N + j other customers in the system. The ra/s satisfy 

m '' = n^ + x-^'"' + ' + xT^'7^ m '- 1 ' im0 - (25) 

The rationale for this equation is that a customer must wait at least 
until the next change of state; the mean of this initial delay is 
(X + Nn)~K If the change of state is caused by an arrival [which 
occurs with probability X/(X + N/i)J, then the customer will have to 
wait an additional period of time whose mean is m J+ i. If, on the other 
hand, the change of state is caused by a departure [which occurs with 
probability Nn/(\ + iV/x)], then with probability j/(j + 1) our 
customer will not be chosen from the group of j + 1 customers, and 
he will have to wait an additional period of time whose mean is wiy-i. 
We now introduce the function 

H(x) = t m-jzi. 

This series converges for x ^ p, since the mean waiting time of an 
arbitrary customer is 

'" = "*W = jvMfeF (26) 
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Multiplying eq. (25) by (1 + p)(j + l)x i and summing on j, we can 
obtain a first-order differential equation : 

*'<*> + ( i- + ,)%"-p) H(x) = h a— )■(»->) • 

The solution to this equation is 

(i _ x \na-f) 2 — x 

7=j) + Wm(1-«)'(8-p) - 

Using eq. (26) as the boundary condition, we see that C must be zero 
in order that H(x) remain finite as x — * p. Thus 

u( \ - 2 ~ x 

H{X) ~ Nfj.il - xY (2 - P )' 

The power series expansion of this function is found to be 
therefore, the means we desire are given by 

i + 2 

vlj Nu(2- P )' 
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